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Multistate  Reliability  Models:  A Survey 

by 

Emad  El-Neweihi  and  Frank  Proschan 

Abstract 

For  a long  time,  the  vast  majority  of  the  models  in  reliability  theory  have 
concentrated  on  the  case  in  which  both  components  and  systems  assume  only  two 
possible  states;  functioning  and  failed.  Unfortunately,  this  represents  a 
gross  oversimplification  of  the  many  real  life  situations  in  which  both  components 
and  systems  actually  assume  a variety  of  states  ranging  from  perfect  operation 
to  complete  failure. 

More  recently,  papers  have  appeared  which  treat  the  more  sophisticated 
and  more  realistic  situations  in  which  components  and  systems  may  assume  many 
states.  In  the  present  paper,  a survey  is  made  of  earlier  work,  but  more 
especially  of  some  quite  recent  work  (some  completed,  some  still  in 
progress)  by  the  relatively  small  number  of  researchers  active  in  this 
important  new  area  of  reliability.  It  is  becoming  apparent  that  this 
research  will  generate  results  not  only  of  value  in  reliability  applications, 
but  also  of  independent  interest  in  multivariate  statistical  analysis. 
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1.  Introduct ion . 


The  vast  majority  of  reliability  analyses  assume  that  components  and 
system  are  in  either  of  the  two  states:  functioning  or  failed.  In  many 
situations  one  is  capable  of  distinguishing  between  various  "levels  of 
performance"  for  both  the  system  and  its  components.  For  such  cases,  the 
existing  dichotomous  model  is  a gross  oversimplification  of  the  real  situation, 
whereas  models  representing  multistate  systems  and  components  are  much  more 
suitable. 

Until  recently,  very  little  work  has  been  done  on  this  more  general 
problem  of  multistate  systems.  Some  earlier  work  treated  only  very  specialized 
aspects  of  multistate  systems,  but  no  comprehensive  treatment  of  these  models 
was  available.  Among  the  earlier  papers  are  [11],  [12],  [13],  [14],  [16], 
and  [17].  However  there  has  been  recently  a growing  interest  in  this  important 
new  area  of  reliability  theory.  More  sophisticated  and  comprehensive  work 
on  multistate  models  has  been  performed  by  Barlow  [2],  El-Neweihi,  Proschan, 
and  Sethuraman  [8]  and  Ross  [15].  In  this  expository  paper  a survey  is  made 
of  the  various  treatments  of  multistate  models.  We  briefly  mention  the 
earlier  work,  but  we  concentrate  on  the  more  recent  and  more  comprehensive 
treatments  of  multistate  models  performed  by  the  relatively  small  number 
of  researchers  active  in  this  important  area  of  reliability. 

We  now  summarize  the  contents  of  this  paper.  Our  terminology  and  notation 
are  similar  to  that  of  Barlow  and  Proschan  [1]  for  the  two  state  case.  In 
Section  2 we  present  the  notation  and  terminology  used  throughout  the  paper. 

In  Section  3 we  consider  a system  of  n components.  For  each  component  and 
for  the  system  itself,  we  can  distinguish  among  different  "levels  of  performance" 
represented  by  a state  space  S.  For  component  i,  x^  denotes  the  corresponding 
state,  i-l,...,n,  the  vector  x = (x^,...,xn)  denotes  the  vector  of  states  of 
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components  l,...,n.  The  state  of  the  system  is  assumed  to  be  a deterministic 
function  of  4>  of  the  states  of  the  components  from  Sn,  the  nth  Cartesian  power 
of  S,  into  S.  Thus  <J>(x,)  is  the  state  of  the  system  corresponding  to  the 
component  state  vector  x-  We  then  survey  the  different  choices  of 
state  space  S,  and  various  definitions  of  the  structure  function  $ presented 
in  different  treatments  of  multistate  models.  We  Investigate  structural 
properties  of  the  various  models,  occasionally  comparing  and  contrasting  them. 

In  Section  4 we  investigate  probabilistic  aspects  of  multistate  models. 

We  survey  the  relationship  (in  a probabilistic  sense)  between  the  performance 
of  the  system  and  the  performances  of  its  components.  For  instance,  system 
performance  is,  as  expected,  a monotone  function  of  component  performances. 

When  the  exact  value  of  system  performance  is  difficult  to  compute,  bounds 
are  provided. 

In  Section  5,  we  survey  dynamic  aspects  of  multistate  systems.  In  earlier 
sections,  it  is  tacitly  assumed  that  time  is  fixed.  In  Section  5,  multistate 
systems  are  viewed  as  operating  over  time.  At  time  0,  the  system  and  each  of 
its  components  are  at  the  maximal  "level  of  performance".  As  time  passes,  the 
performance  levels  of  components  (and  consequently  of  the  system)  deteriorate 
to  successively  lower  levels  until  finally  level  0 (complete  failure)  is 
reached.  Concepts  of  IFRA  and  NBU  stochastic  processes,  analogous  to  the 
corresponding  lifelength  distributions  in  the  binary  case  are  defined  and 
studied  by  various  researchers . Some  generalized  IFRA  and  NBU  closure  theorems 
are  presented. 

Finally  in  Section  6,  we  show  by  means  of  two  examples  how  theories  of 
multistate  systems  may  be  applied  to  existing  binary  reliability  models. 


r 


2.  Notation  and  Terminology. 

The  vector  x = (x^,...,xn)  denotes  the  vector  of  states  of  components 

1) • • « jU* 

C - {l,2,...,n}  denotes  the  set  of  component  indices. 

— (x^, . • • * • • »J  * * * **n)  w^ere 

(*i,x)  = 0^,...,  x±_v  *»  x1+1»***»xn> 

j_= 

■<j_  < x means  that  y^  £ x^  for  i-1 n and  y^  < x^  for  some  i 

a * is  a probability  vector  means  that  0 £ £ 1, 


M 

j-0,1, . . . ,M  and  Z a 


j-o 


j 


a £ a' ' , where  both  a,  a'  are  probability  vectors,  means  that 
st 


Z ci^  £ Z a^,  j-0,1,. 
k-j  K k-j  K 


,M. 


xVy  denotes  max(x,y). 
xVy.  = (x1Vy1 .... »xnVyn) . 
xAy  denotes  min(x,y). 
xAy_  = (x1Ay1, . . . • 

"Increasing"  is  used  in  place  of  "nondecreasing"  and  "decreasing" 

is  used  in  place  of  "nonincreasing".  When  we  say  f(x^,...,xQ) 
is  increasing  we  mean  f is  increasing  in  each  argument. 

Given  a univariate  distribution  F,  its  complement  1-F  is  denoted  by  F. 
Given  a set  S,  Sn  denotes  its  nth  Cartesian  power. 
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3.  Deterministic  Properties  of  Multistate  Coherent  Systems. 

Consider  a system  of  n components.  We  assume  that  the  performance  of  the 
system  depends  deterministically  on  the  performances  of  the  components.  Thus 
given  x,  the  vector  of  component  states,  we  may  determine  $ (x) , the  system 
state.  The  function  <J>  is  called  the  structure  function  of  the  system.  In 
the  binary  case,  it  is  assumed  that  both  components  and  system  are  in  either 
of  two  states:  functioning  or  failed.  The  variables  x^,  i-l,...,n,  as  well 
as  <J>(x),  assume  their  values  in  the  state  space  S ■ {0,1},  where  0 denotes 
failure  and  1 denotes  functioning.  The  structure  function  <J>  is  then  a map 
from  {0,l}n  into  {0 , 1 } . The  structure  function  <p  satisfies  certain  conditions 
that  represent  intuitively  reasonable  properties  of  systems  encountered 
in  practice.  The  following  two  conditions  are  required  for  a binary  system  to 
be  a coherent  structure  ([1],  Def.  2.1,  p.  6), 

(i)  The  function  c|>(x)  is  increasing. 

(ii)  Each  component  is  relevant  to  the  system,  i.e. , for  each 
i there  exists  a vector  (*i,x)  such  that  <j>(l^,jc)  > <j>(0^,x).  This 
means  that  the  function  <{>  is  not  constant  in  its  ith  argument, 

1st  1 , . . . , n . 

Condition  (i)  embodies  the  reasonable  assumption  that  improving  the  performance 
of  a component  is  not  harmful  to  system  performance.  Condition  (ii)  eliminates 
from  consideration  components  which  have  no  effect  on  system  performance.  The 
theory  of  binary  coherent  structures  has  served  as  a unifying  foundation  for 
a mathematical  and  statistical  theory  of  reliability  for  the  dichotomous  case. 

The  binary  model,  however,  is  an  oversimplification  in  describing  a 
situation  in  which  either  the  system  or  its  components  (or  both)  are  capable 
of  assuming  a whole  range  of  levels  of  performance,  varying  from  perfect 
functioning  to  complete  failure.  In  these  situations,  models  representing 
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mulitstate  systems  and  multistate  components  are  much  more  useful  in  describing 
system  performance  in  terms  of  component  performances.  Naturally,  the  first 
step  in  constructing  such  models  is  to  provide  useful  definitions  of  state 
spaces,  representing  the  sets  of  levels  of  performance,  and  of  the  structure 
function  <}>,  that  relates  the  performance  of  the  system  to  the  performances  of 
its  components.  A theory  of  multistate  structures  can  then  serve  as  a unifying 
foundation  for  a mathematical  and  statistical  theory  of  reliability  in  the 
multistate  case.  Among  the  earlier  attempts  to  this  type,  we  mention  the 
following  two  examples: 

Hirsch,  et  al  [11],  in  a treatment  of  "cannibalization",  consider  a system 
of  n components;  the  state  of  component  i is  represented  by  the  binary 
variable  x^  assuming  the  values  1 or  0 according  to  whether  component  i is 
functioning  or  failed,  1*1,..., n.  However  the  system  itself  can  be  in 
any  of  M4-1  states  representing  various  levels  of  performance.  The  set  of 
possible  performance  levels  is  assumed  to  be  totally  ordered  and  is  then 
represented,  without  loss  of  generality,  by  S * {0,1, ...,M},  where  0 denotes 
complete  failure  and  M denotes  perfect  performance.  The  structure  function 
<{>  is  a map  from  {0,1}  into  S.  Note  that  S can  have  at  most  2 elements. 

The  structure  function  4>  is  assumed  to  be  monotone  increasing,  with  <})(0)  * 0, 
4>(1)  * M.  The  authors,  however,  do  not  attempt  a general  treatment  of 
multistate  models.  Their  main  concern  rather  is  to  investigate  the  mathematical 
model  for  cannibalization  to  determine  how  components  may  be  exchanged  to 
improve  system  performance.  The  following  multistate  model  is  presented  by 


Postelnica  [14]:  Consider  a system  of  n components  in  which  the  state  space 
for  the  system  and  for  each  of  its  components  is  the  unit  interval  [0,1], 


representing  a continuous  range  of  performance  from  perfect  performance  (1) 
to  complete  failure  (0).  The  structure  function  <f>:  [0,l]n  -*■  [0,1]  satisfies 
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the  following  conditions: 


(a) 

MI)  - l. 

(b) 

<J>(0)  = o. 

(c) 

<]>(x)  is  monotone  increasing. 

(d) 

Me)  >.  c. 

0 < c < 1. 

(e) 

Mc^,0)  < c 

0 < c < 1. 

The  author  does  not  attempt  a comprehensive  treatment  of  such  multistate 

structures,  buth  rather  investigates  some  very  special  applications. 

More  recent  and  more  comprehensive  research  in  multistate  systems  has 

been  performed  by  Barlow  [2],  El-Neweihi,  Proschan  and  Sethuraman  [8]  (hereafter 

referred  to  as  EPS  [8]),  and  Ross  [15].  The  definition  given  by  Barlow  [2] 

for  the  multistate  structure  is  set-theoretical,  based  on  the  concept  of  min 

path  sets  and  min  cut  sets  of  binary  coherent  structures.  Consider  a system 

of  n components.  Assume  that  the  state  space  for  each  of  the  components  as 

well  as  for  the  system  is  the  set  S = (0,1, ...,M},  where  0 denotes  the  failed 

state  and  M denotes  the  maximal  or  perfect  state.  Let  P^,...,Pr  be  non-empty 

r ■ 

subsets  of  C such  that  U P.  = C and  P.  P . , i^j . The  structure  function 

i=l  1 J 

<p:  Sn  S is  defined  by 


<j>(x) 


max  min  x^ , 


(3.1) 


l£j<r  iePj 

where  xeSn  is  the  vector  representing  the  states  of  components  l,2,...,n. 

In  the  binary  case  the  structure  function  given  in  (3.1)  is  the  most  general 
coherent  structure  ([1],  Chapter  1),  and  the  sets  P^,...,Pr  are  its  min  path 
sets.  Let  <J>'  be  the  binary  coherent  structure  associated  with  P^,...,Pr<  The 
multistate  coherent  structure  <J>  specified  in  (3.1)  can  then  be  expressed  in 
terms  of  the  corresponding 


each  1*1,..., n.  let  y 


U 


binary  coherent  structure  <J>'  as  follows:  For 
' 1 if  x±>j 

" '“U 'nj 


1° 


o .w. 


and  let  “ (y-n  » • • • *y„* ) » j“0,l, . . . ,M. 
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It  is  fairly  easy  to  see  that  4>(x)  >.  j iff  <j>'(£.)  * 1,  and 

M 

<p(x)  - E 4>'(Z .).  (3.2) 

J-l  3 

Thus  the  multistate  coherent  structure  given  by  Barlow  [2],  is  very 
closely  related  to  a corresponding  binary  coherent  structure.  Exploiting 
this  relationship  makes  it  easy  to  extend  results  from  the  binary  case  to  the 
multistate  case. 

A more  general  approach  has  been  taken  by  EPS  [8]  to  define  multistate 
coherent  structures.  The  common  state  space  for  each  of  the  components  and 
for  the  system  is  again  taken  to  be  the  set  S = {0,...,M},  representing  the 
Mfl  levels  of  performance  ranging  from  complete  failure  (0)  to  perfect 
functioning  (M) . The  structure  function  <j>:  Sn  -*■  S is  assumed  to  satisfy  three 
conditions. 

Definition  3.1.  A system  of  n components  is  said  to  be  a multistate 
coherent  system  (MCS)  if  its  structure  function  <}>  satisfies: 

(i) '  <{>  is  increasing. 

(ii) '  For  level  j and  component  i,  there  exists  a vector 
(*i,x)  such  that  (Kj^x)  = j while  ^(i^.x)  + j for  j , i-l,...,n, 
and  j=0,l,. . . ,M. 

(iii)’  <K±)  = j for  j=0,l,...,M. 

The  three  axioms  embodied  in  Definition  3.1  extend  the  notion  of 
binary  coherent  system  to  the  new  notion  of  a multistate  coherent  system. 

Note  that  conditions  (i) ' and  (ii)'  generalize  conditions  (i)  and  (ii) 
in  the  binary  case.  Condition  (iii)'  is  automatically  satisfied  in  the 
binary  case,  but  is  not  implied  in  the  present  multistate  case  by  (i) ' 
and  (ii)'.  Also  note  that  since  the  structure  function  in  (3.1),  defined 
by  Barlow  [2],  satisfy  conditions  (i) ' , (ii)',  and  (iii)'  of  Definition 
3.1,  they  constitute  a subclass  of  the  MCS  class.  For  instance,  for  a two 


component  system,  only  two  distinct  systems  satisfy  (3.1),  namely,  the  parallel 
and  the  series  system,  regardless  of  the  cardinality  of  S.  However  for 
S = {0,1,2},  there  are  more  than  12  MCS's. 

The  definition  given  by  Ross  [16]  for  a multistate  system  is  less  structured 
than  either  the  Barlow  [2]  specification  or  the  EPS  [8]  specification.  The 
state  space  S is  taken  to  be  [0,°°)  and  the  structure  function  <}>  is  any  monotone 
increasing  function  from  [0,°°)n  into  [0,°°).  Ross  [15]  has  not  attempted  to 
investigate  structural  properties  of  his  model;  rather,  he  concentrates  on  the 
stochastic  properties  of  his  model  when  observed  either  at  a fixed  point  in 
time,  or  when  observed  at  different  points  in  time  (dynamic  models).  Results 
of  this  type  will  be  surveyed  in  the  next  two  sections. 

In  the  remainder  of  this  section  we  present  various  structural  properties 
of  the  multistate  structures  given  by  Barlow  [2]  and  EPS  [8].  These  properties 
extend  well-known  results  in  the  binary  case  ([1],  Chapter  1)  to  the  more 
general  multistate  case. 

The  following  theorem  gives  simple  bounds  on  MCS  performance: 

Theorem  3.1.  Let  4>  be  the  structure  function  of  an  MCS  of  n components. 

Then 

min  x.  £ <J>(x)  £ max  x..  (3.3) 

i£i£n  l£i£n 

Theorem  3.1  states  that  a parallel  system  yields  the  best  performance 
of  an  MCS,  and  a series  system  yields  the  worst  performance.  Using  this 
theorem,  EPS  [8]  establish  probabilistic  bounds  on  system  reliability. 

As  in  the  binary  case,  the  following  lemma  in  EPS  [8]  gives  a decomposition 
identity  useful  in  carrying  out  inductive  proofs.  It  holds  for  any  multistate 
structure,  not  just  for  the  MCS. 
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Lemma  3.1.  The  following  identity  holds  for  any  n-component  structure 
function  <J>: 

M 

<}>(x)  - E for  (3.4) 

j-0  lxi  31 

where 

1 if  xi=j 

[*i*j  J o if  # 

As  in  the  binary  case,  EPS  [8],  define  a dual  structure  for  each  multistate 
structure. 

Definition  3.2.  Let  <p  be  the  structure  function  of  a multistate  system. 

The  dual  structure  function  <j)D  is  given  by: 

<t>D(x)  =*  M - (|) (M-x^ , M-x2,...,  M-xn>.  (3.5) 

It  is  easy  to  verify  that  the  dual  of  an  MCS  is  an  MCS. 

Example  3.1.  The  dual  of  a series  (parallel)  system  is  a parallel  (series) 
system.  More  generally,  the  dual  of  a k-out-of-n  is  an  (n-k+l)-out-of-n 

system,  where  a k-out-of-n  system  is  given  by  <{>(x)  = x(n_jc+i)* 

Design  engineers  have  used  the  well-known  principle  that  redundancy  at 
the  component  level  is  preferable  to  redundancy  at  the  system  level.  This 

principle  is  presented  by  EPS  [8]  in  mathematical  form  in  (i)  of  the 

following  theorem;  (ii)  is  a dual  result. 

Theorem  3.2.  Let  | be  a structure  function  of  an  MCS.  Then 

(i)  <p  (xVjO  _>  <p(x)V<t>  (jr)  . 

(ii ) <p  (xAy)  1 <J>  (x) A<t>  (z)  • 

Equality  holds  in  (i)  for  all  x an^  Z ff  an<^  only  if  the  structure  is  parallel. 
Equality  holds  in  (ii)  for  all  x and  if  and  only  if  the  structure  is  series. 
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Parts  of  (i)  and  (ii)  of  Theorem  3.2  are  also  proved  by  Barlow  [2]. 

In  the  binary  coherent  structures  the  concepts  of  minimal  path  vectors 
and  minimal  cut  vectors  play  a crucial  role.  The  analogue  in  MCS  theory  is 
the  concept  of  critical  connection  vectors.  This  concept  is  defined  by 
EPS  fS]  in  the  following: 

Definition  3.3.  A vector  x is  said  to  be  a connection  vector  to  level 
I if  <}>(x)  =*  j,  j*0,l, . . . ,M. 

Definition  3.4.  A vector  x is  said  to  be  an  upper  critical  connection 
vector  to  level  j if  <J>(x)  * j and  £ < x implies  <K2.)  < j»  J"1»...»M. 

A lower  critical  connection  vector  to  level  j can  be  defined  in  a dual 
manner,  j=0,l, . . . ,M-1. 

The  existence  of  such  critical  connection  vectors  is  guaranteed  by  the 
conditions  of  Definition  3.1. 

Let  x be  an  upper  critical  connection  vector  to  level  j . Define 
C j (x)  = {i:  xi>(j}.  Obviously  Cj(x)  is  a non-empty  subset  of  C * {l,...,n}. 

For  j*»l,...,M,  let  = {C^(x):  x is  an  uPPer  critical  connection  vector  to 
level  j}.  Then  the  following  lemma  by  EPS  [8],  shows  that  enjoys  a 

property  similar  to  that  enjoyed  by  the  minimal  path  sets  and  minimal  cut  sets 
in  the  binary  case. 


Lemma  3.2.  For  j * 1 , . . . ,M, 


UCj  = {1, . . . ,n}. 

For  j=l,...,M.  let  , . . . he  the  upper  critical  connection  vectors 

J 

to  level  j,  where  > • • • »y„r^  ’ 1— r— n j ‘ 1116  followin8  theorem  by  EPS  [8], 

expresses  the  state  of  an  MCS  using  its  upper  critical  connection  vectors. 

Theorem  3.3.  Let  <f>  be  the  structure  function  of  an  MCS.  Let 
be  its  upper  critical  connection  vectors  to  level  j,  j*l,...,M.  Then 
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<Kx)  j if  and  only  if  x > for  some  j £t  <_  M and  some  l£&<nt. 

The  above  theorem  Is  uitllized  to  establish  bounds  on  the  system  performance 
distribution,  as  will  be  shown  In  the  next  section. 

4.  Probabilistic  Properties  of  Multlstate  Coherent  Systems. 

The  deterministic  relationships  between  the  performance  of  a multistate 
system  and  that  of  Its  components  are  exploited  by  the  various  researchers  in 
the  field  to  investigate  the  probabilistic  properties  of  multistate  systems. 

In  this  section  we  survey  important  relationships  between  the  stochastic 
performance  of  the  system  and  the  stochastic  performances  of  its  components. 

These  results  provide  bounds  on  system  performance  which  are  particularly 
useful  when  exact  system  performance  is  difficult  to  evaluate. 

Let  X^  denote  the  random  state  of  component  i,  i«l,...,n.  Let  X - (X^,...,Xn) 
be  the  random  vector  representing  the  states  of  components  l,...,n,  where  the 
X^'s  are  assumed  to  be  stochastically  mutually  independent.  Then  <J>(X)  is  the 
random  variable  representing  the  system  state,  where  <J>  is  the  structure  function 
of  the  system.  Naturally,  the  random  variables  assume  their  values  in  the  state 
space  S according  to  certain  probability  laws.  In  the  model  described  by 
Postelnicu  [14],  X^,...,Xn  as.  well  as  <J>(X)  are  distributed  in  the  unit  interval 
[0,1],  with  cumulative  distribution  functions  F^,...,Fn  and  F respectively. 
Postelnicu  [14]  discusses  briefly  bounds  on  F in  terms  of  F^,...,Fn.  In  the 
models  described  by  Barlow  [2]  and  EPS  [8],  the  random  variables  X^,...,Xn 
and  <j>(X)  assume  their  values  in  S * {0,...,M},  with 


j-0,1. 


P[X±-j ] - Pij 
P[X±<J ] - Pi(j) 
P[<KX)-J]  - Pj 
P[«J>(X)<J]  - P(j), 


(4.1) 


,M,  and  i-l,2,...,n.  P^  represents  the  performance  distribution  of 
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component  i,  while  F represents  the  performance  distribution  of  the  system. 
Clearly, 

j 

Mj)  - ^ p 
k*0  X 

Pi(M)  = 1, 

for  1*1,..., n.  Similar  relationships  hold  for  P.  Let  h * E<|i(X);  we  may 
express  h as  follows: 

h — h (P.»...,P  ), 

£ -i  n 

since  h is  a function  of  the  P^,...,Pn.  Alternatively,  we  may  epxress  h as 


follows : 


h ” hp(£.^»  • • • i 


where  jj  = (p^g,  p^,...,  piM)  for  i*l,...,n.  In  either  case,  EPS  [8]  calls 
h the  performance  function  of  the  system. 

Using  Lemma  3.1,  EPS  [8],  expresses  the  performance  function  of  a system 
of  components  in  terms  of  performance  functions  of  systems  of  n-1  components. 
Such  a decomposition  identity  is  useful  in  carrying  out  a proof  by  induction 
and  in  deriving  properties  of  h. 

Lemma  4.1.  The  following  identity  holds  for  h: 

M 

h (p^> « * » >Pq)  “ ^ P^j  h(j^,  P^ » • • • »P^) » i“l» • • • »n»  (4.2) 

where  h(Ji,  £1»**.,2n)  = E 

The  following  theorem  due  to  EPS  [8]  shows  that  h is  strictly  increasing 
in  each  p^ , j>0.  This  property  generalizes  the  corresponding  well  known 
property  of  h in  the  binary  case. 


! 


T 








— 
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Theorem  4.1.  Let  h^,...,^)  be  the  performance  function  of  an  MCS. 

Let  0<pi^<l  for  1-1,..., n,  j-0,1, . . . ,M.  Then  is  strictly 

increasing  in  p^  , i"l n,  j-l,...,M. 

Properties  of  h as  a function  of  P^.-.jP^  are  also  investigated  by 
Barlow  [2]  and  EPS  [8].  The  following  theorem  due  to  EPS  [8],  shows  that 
h(P^,...,Pn)  is  monotone  increasing  with  respect  to  stochastic  ordering. 

A similar  result  is  proved  by  Barlow  [2]  for  his  subclass  of  the  MCS  (see 
(3.1))  using  a different  proof.  The  same  property  is  also  proved  by  Ross 
[15]  for  his  multistate  model. 

Theorem  4.2.  Let  P^,  P^  be  two  possible  performance  distributions 
for  component  i,  1-1,. ..,n.  Assume  P^(j)  >^P^(j)  for  j“0,l,...,M, 
i-l,...,n.  Let  P(P')  be  the  corresponding  system  performance  distribution. 
Then 

(i)  P(j)>P'(j)  for  j-0,1,... ,M, 

(ii)  h(P. P)  < h(P’ P’).  (4.3) 

i n — -L  n 

A useful  decomposition  identity  is  given  by  EPS  [8]  for  P[<j>(x)>£]  , 
namely 

Theorem  4.3.  Let  4>  be  a multistate  structure  function.  Then 

M 

P[<J>(X)>A]  - Z p,.  P[<Mj  ,X)>£],  JL=1,...,M.  (4.4) 

j-0  13 

Relation  (4.4)  expresses  the  survival  probability  of  a structure  of 
n components  in  terms  of  survival  probabilities  of  structures  of  n-1 
components . 

Using  Theorem  3.1,  EPS  [8]  obtain  the  following  useful  bounds  on  P 

and  h in  terms  of  P, ,...,P  : 

1 n 

Let  P be  the  performance  distribution  andh  be  the  performance  function 
of  an  MCS.  Let  P^  be  the  ith  component  performance  distribution,  i-l,...,n. 
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Then  for  j-0,1, . . . ,M-1: 


n n 

IT  P,(j>  1 P(j)  £ 1 - Tt  P,  0), 
i-1  1 i-1 

H n M n 

Z ir  P.  (j-1)  £ h < Z [1-  it  P.  ( j— 1] , 

j-1  i-1  1 j-1  i-1  1 

where  P^Cj)  - 1 - Pj^j). 

The  concept  of  upper  connection  critical  vectors  introduced  by  EPS  [8] 
is  exploited  to  establish  further  bounds  on  P and  h.  Let  be  ttie 

upper  critical  connection  vectors  to  level  j,  j-l,...,M  (see  Definition 

3.4) .  Let  denote  the  event  [X  :>  r-l,...,^.  By  Theorem  3.3, 

M °t 

P[<KX)>j]  - P[  U U A ].  Now  using  the  well  known  inclusion-exclusion 
t*j  r-1  r 

principle,  the  authors  establish  upper  and  lower  bounds  on  P[<J>(X)>j]  - P(j-l). 
Note  that  P(A^.)  - PTX  >.  ptx^  — for  ^ r— nj  and  j_1»"*»M* 

An  interesting  generalization  of  the  Moore-Shannon  Theorem  [[1],  Theorem 

5.4]  is  obtained  by  Barlow  [2]  using  the  close  relationship  between  his 
definition  of  a multistate  coherent  system  and  that  of  the  binary  coherent 
system.  Recall  that  corresponding  to  every  multistate  structure  function  <{> 
defined  by  Barlow  [2],  there  is  a binary  coherent  structure  4> * closely 
related  to  <J>  (see  (3.1)  and  (3.2)).  Let  h'  be  the  binary  reliability  function 
associated  with  <J>' , i.e.,  h'  - E<J>'(Y),  where  Y - (Y^,...,Yn)  is  a random 
vector  whose  components  are  Bernoulli  random  variables.  In  view  of  (3.2), 

it  is  easily  verified  that 

(4.5) 


P[4>(X)>j]  = E^’Ofj)  - h'  (cLj), 

M 

where  ^ - (q^  , . . . ) , and  * Z Pife,  i=l,...,n. 

k— j 

Recall  that  Moore  and  Shannon  show  that  binary  coherent  reliability 
functions  are  S-shaped  in  the  sense  that  if  all  components  function  with 
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probability  p,  either  h(p)  > p or  h(p)  £ p for  all  0<p£l,  or  there  exists 
0<pQ<l  such  that  h(p)  £ p for  0<p£pQ,  while  h(p)  >_  p for  l^p>pQ-  Barlow  [2] 
gives  a natural  generalization  of  this  result  to  the  multistate  case  with 
respect  to  stochastic  ordering. 

Theorem  4.4.  Let  “ * (<*0 0^)  for  1*1,..., n.  Assume  h'(pQ)  - pQ 

* 

(0<pQ<l).  Let  a * (l-pQ,0,. . . ,0,Pq).  Then 

* 

(a)  a £ a implies  that  £ £ a, 

st  st 

* 

(b)  a £ a implies  that  £ >^  ct, 

st  st 

where  p - (p  , ...,p„),  P.  - p[<f>(X)-j],  j*0,...,M,  and  a'  < a"  means  that 
M ~ M ° M J st 

E CL*  £ E CL  , j*0 ,1, . . . ,M. 
k*j  k-3 

Note  that  (4.5)  is  central  to  the  proof  of  the  above  theorem. 

Finally,  in  the  model  proposed  by  Ross  [16],  X^,  i*l,...,n,  and 

<J>(X)  are  non-negative  random  variables  with  distribution  functions  F^, 

1*1,..., n,  and  F respectively.  The  function  r(F^,...,Fn>  is  defined  by 

r(F1,...,F  ) - E4»(X), 
in  — 

where  F^  * 1-F^,  •*tn* 

Using  an  extension  of  Lemma  2.3,  p.  84,  of  Barlow  and  Proschan  [1], 

Ross  [15]  proves  the  following: 

Theorem  4.5.  If  <{»  is  a binary  function  then 

r(fj,...,fj[)  > [r^ Fn)]a  (4.6) 

for  all  0£a£l. 

As  a consequence  of  che  above  theorem,  Ross  [15]  proves: 

Corollary  4.1.  Let  X^,...,X  be  independent  IFRA  random  variables. 


Then 
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n 

(a)  I X,  is  IFRA. 
i-1 

n n 

(b)  P{  7T  X.  > aa  } > (P{  7T  X.  > a})  , 0<a<l. 

i-1  i-1 

Recall  that  a distribution  function  F with  F(0)  - 0 is  said  to  be  an 
increasing  failure  rate  average  (IFRA)  distribution  if 

F(otx)  ^ [F(x)]a  for  all  0<tt<l,  x>0. 

Observe  that  part  (a)  of  Corollary  4.1  represents  the  well-known  property 
of  the  closure  of  the  IFRA  distributions  under  the  convolution  operation. 

Ross  [15]  also  utilizes  Theorem  4.5  in  proving  a generalized  IFRA 
closure  theorem  which  is  presented  in  the  next  section. 

5.  Dynamic  Models  for  Multistate  Coherent  Systems. 

In  previous  sections,  we  consider  deterministic  and  probabilistic 
properties  of  multistate  systems  at  a fixed  point  in  time.  In  this  section 
we  survey  some  dynamic  aspects  of  multistate  structures.  We  now  consider 
multistate  system  as  operating  over  time.  At  time  0 the  system  and  each  of 
its  components  are  at  their  maximal  level  of  performance.  As  time  passes, 
the  performance  levels  of  components  (and  consequently  of  the  system) 
deteriorate  to  lower  levels  until  finally  level  0 (complete  failure)  is 
reached . 

In  the  binary  case,  the  length  of  time  during  which  a component  (system) 
functions  is  called  the  llfelength  of  the  component  (system);  each  lifelength 
is  a non-negative  random  variable.  The  corresponding  lifelength  distribution 
has  been  classified  according  to  various  notions  of  aging.  See,  e.g. , [1]. 

Two  of  the  important  classes  of  life  distributions  are  the  increasing  failure 
rate  average  (IFRA)  class  and  the  new  better  than  used  (NBU)  class.  Closure 
of  these  classes  under  various  basic  reliability  operations,  such  as 
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convolution  of  distributions  and  formation  of  binary  coherent  systems,  is 
demonstrated  in  [1],  The  counterparts  of  these  concepts  in  the  multistate 
case  have  been  Investigated  by  Barlow  [2],  EPS  [8],  and  Ross  [15]. 

Let  (X^(t),  t>0}  denote  the  decreasing  stochastic  process  representing 
the  state  of  component  i at  time  t,  where  t ranges  over  the  non-negative 
real  numbers  for  i»l,...,n.  The  stochastic  process  {<J>(X(t)),  t>0}  is  also 
decreasing  and  represents  the  corresponding  system  state  as  time  varies, 
where  X(t)  * (X^(t) , . . . ,Xn(t) ) . The  processes  (X^(t),  t>0},  i=l,...,n  are 
assumed  to  be  mutually  independent. 

In  Barlow's  model  where  the  state  space  is  {0,1,...,M},  let  us  call 
{ j , j+l,...,M}  the  "good"  states.  Assume  that  [P{x^(t)  is 

decreasing  in  t>0  for  fixed  j.  It  is  easily  verified  that  [P{<{>(X(t))  j}]^*" 

is  decreasing  in  t>0  for  fixed  j . Thus  the  above  result  states  that  if  the 
length  of  time  spent  by  each  component  in  the  "good"  states  is  an  IFRA 
random  variable,  then  the  corresponding  length  of  time  spent  by  the  multi- 
state system  in  the  "good"  states  is  also  an  IFRA  random  variable.  In  the 
binary  case  this  represents  the  so-called  IFRA  closure  (under  formation 
of  binary  coherent  systems)  theorem.  Note  that  from  (4.5)  the  proof  of 
the  IFRA  closure  theorem  for  Barlow's  model  is  immediate. 

The  following  definition  is  due  to  Ross  [15]. 

Definition  5.1.  The  stochastic  process  {X(t),  t>0}  is  said  to  be  an 
IFRA  process  if  T&  * inf{t:  X(t)  _<  a}  is  an  IFRA  random  variable  for 
every  a ^ 0. 

Having  introduced  this  definition,  Ross  [15]  then  proves  the  following 
generalized  IFRA  closure  theorem. 

Theorem  5.1.  Let  (X^t),  t>0},  i»l,...,n,  be  independent  IFRA 
processes  and  <j>  a multistate  structure  function.  Then  {<j>(X(t)),  t>0}  is 
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an  IFRA  process. 

The  crucial  tool  in  proving  the  above  theorem  is  Theorem  4.5. 

Ross  [15]  also  defines  an  NBU  process  and  proves  a generalized  NBU 
closure  theorem.  First  let  us  recall  the  definition  of  an  NBU  random 
variable. 

Definition  5.2.  A non-negative  random  variable  Y with  distribution 
function  F is  said  to  be  new  better  than  used  (NBU)  if  F(s+t)  F(s)F(t) 
for  all  s>0,  t>0,  where  F - 1-F. 

Now,  Ross  [15]  gives  the  following  definition  of  an  NBU  process. 

Definition  5.3.  The  decreasing  stochastic  process  {X(t)  t>0}  is  said 
to  be  NBU  if  with  probability  1, 

P{Ta>s+t|X(u)>a,  0 <u<s } < P(Ta>t} 
for  all  s,t,  a>0,  where  T is  as  in  Definition  5.1. 

“ a 

Using  his  definition,  Ross  proves: 

Theorem  5.2.  If  the  component  processes  are  independent  NBU  processes, 
then  {<f>(X(t)>,  t>0)  is  also  NBU. 

Another  definition  of  an  NBU  process  is  given  by  EPS  [8],  and  then  a 
simple  characterization  for  this  NBU  process  is  derived.  Using  their 
characterization,  they  give  a simple  proof  of  a generalized  NBU  closure 
theorem.  The  EPS  definition  of  an  NBU  process  .is  as  follows: 

Definition  5.4.  The  stochastic  process  {X^Ct),  t>0)  is  an  NBU  stochastic 
process  if  T,  . = inf{t:  X. (t)  < j } is  an  NBU  random  variable  for  j-0,...,M-l. 

Recall  that  the  state  space  for  the  EPS  [8]  model  is  the  set  {0,...,M}. 

The  following  lemma  gives  a simple  characterization  of  an  NBU  process. 

Lemma  5.1.  The  stochastic  process  (X^t),  t>0)  is  NBU  if  and  only  if 
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st 

Xi(8+t)  < min(X^(s),  X'(t», 

where  X^(s)  and  X^(t)  are  two  Independent  random  variables  having  the  same 
distributions  as  X^(s),  X^(t)  respectively. 

Using  their  Lemma  5.1,  EPS  [8],  prove  the  following  generalized  NBU 
closure  theorem. 

Theorem  5.3.  Let  <j>  be  the  structure  function  of  an  MCS  having  n 
components  and  (X^(t),  t>0}  be  the  ith  component  performance  process, 
l>l,...,n.  Let  {X^Ct),  t>0},  i*l,...,n,  be  mutually  independent  NBU  processes. 
Then  {<£(X(t)),  t>0}  is  an  NBU  stochastic  process. 

Remark  5.1.  The  useful  characterization  of  Lemma  5.1,  adapted  to  the 
binary  case,  yields  a simpler  proof  of  the  NBU  closure  theorem  than 
the  proof  given  in  [1]. 

6.  Applications  of  Multistate  Reliability  Models. 

In  this  section  we  illustrate  by  means  of  two  examples  that  the  theory 
of  multistate  reliability  models  provides  useful  and  new  treatments  of  some 
existing  binary  reliability  models.  This  shows  that  not  only  do  the  multi- 
state reliability  models  provide  more  realistic  analyses  of  many  real  life 
situations,  but  they  also  permit  us  to  obtain  a better  understanding  and 
a more  efficient  treatment  of  existing  models  in  the  two-state  case. 

Example  1 appears  in  El-Neweihl , [ 9 ] . 

Example  1;  EPS  [7]  study  the  following  model.  A series-parallel  system 
consists  of  k+1  subsystems  Cq,C^, . . . ,0^,  also  called  cut  sets.  Cut  set 
contains  n^  components  arranged  in  parallel,  i«0,l,...,k.  No  two  cut  sets 
have  a component  in  common.  Components  fail  one  at  a time,  and  after  t 
components  have  failed,  each  of  the  remaining  components  is  equally  likely 
to  fail,  t*0,l,...  . The  system  fails  upon  failure  of  any  of  the  cut  sets; 
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cut  set  fails  when  all  of  its  components  fail.  This  model  has  many  applications 
in  the  study  of  reliability,  extinction  of  species,  inventory  depletion,  urn 
sampling,  among  others. 

In  part  I of  [7],  EPS  study  the  probability  P(ng;n)  that  the  system  fails 
because  of  a specified  cut  set,  say  Cq,  fails  first.  Several  alternative 
expressions  and  recurrence  relations  for  this  probability  are  obtained.  Some 
of  these  formulae  are  useful  in  the  computation  of  desired  quantities,  while 
others  are  used  to  demonstrate  qualitative  features  like  monotonicity,  Schur- 
concavity,  etc.,  and  derive  asymptotic  limits.  Similar  results  are  also 
obtained  for  the  more  general  model  in  which  an  "alarm"  rings  when  a cuts  set 
size  first  reaches  a,  where  a is  a specified  positive  integer. 

In  part  II  of  [7],  the  authors  study  the  probability  distribution, 
frequency  function,  and  failure  rate  of  the  lifelength  of  series-parallel 
systems,  where  system  lifelength  refers  to  the  number  of  components  that 
have  failed  at  the  time  of  the  system  failure. 

We  now  show  the  relationship  between  the  above  model  and  the  multistate 

models.  Assume  n^=*M,  i=0,...,k.  Let  the  state  of  component  i of  an  MCS 

be  defined  to  be  the  number  of  functioning  components  in  cut  set  C^.  Thus 

the  above  model  may  be  viewed  as  a series  MCS.  Let  X^,...,X^  be  the  random 

variables  representing  the  lifelengths  of  components  in  cut  set  C^,  i«0,l,...,k. 

For- multistate  component  i,  the  lengths  of  time  spent  in  state  M,  M-l,...,l 

are  given  by  X*^,  X^  - X*  ,...,Xq^  - X*M_^,  i=0,...,k,  where 

X1  X1  X.(1)  i i 

A(l)’  a(2)  * * * * are  the  M order  statistics  of  X^,...,X„. 

Such  an  identification  relating  the  multistate  model  and  the  binary  model 
permits  us  to  answer  a host  of  questions  concerning  one  of  the  two  models 
using  results  obtained  for  the  other.  For  instance  one  can  find  the  probability 
that  component  0,  say,  reaches  an  "alarm"  state  j,  say,  first.  Such  information 
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Is  helpful  In  planning  maintenance  and  replacment  policies 

Example  2.  Consider  a series  system  of  n binary  components.  Assume  we 
have  M-l  spares  for  each  of  the  n components.  A failed  component  is 
instantaneously  replaced  by  one  of  its  spares.  When  the  original  component 
1 is  functioning  (and  thus  none  of  the  spares  has  been  used),  we  consider 
that  component  i is  in  state  M.  Upon  failure  of  the  original  component 
one  of  its  spares  is  used  to  replace  it,  and  so  the  component  now  enters 
state  M-l,  etc.  Thus  we  can  view  the  system  with  its  spares  as  a multi- 
state series  system.  Let  X^,...,X^be  random  variables  representing  the 
lifelengths  of  component  i and  of  its  spares,  i=l,...,n.  Assume  that  all 
the  random  variables  are  mutually  independent.  Obviously,  the  length  of 
time  spent  by  component  i in  a single  "state"  or  in  a group  of  "states" 

can  be  expressed  in  terms  of  a sum  of  an  appropriate  subset  of  X^ X^, 

1=1,..., n.  Thus  we  may  view  a binary  system  with  spares  as  a multistate 
system.  Again,  such  an  identification  is  mutually  beneficial  in  the  study 
of  both  models. 
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states.  In  the  present  paper,  a survey  is  made  of  earlier  work,  but  more 
especially  of  some  quite  recent  work  (some  completed,  some  still  in 
progress)  by  the  relatively  small  number  of  researchers  active  in  this 
important  new  area  of  reliability.  It  is  becoming  apparent  that  this 
research  will  generate  results  not  only  of  value  in  reliability 
applications,  but  also  of  independent  interest  in  multivariate  statistical 
analysis. 
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